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Abstract. We describe and model emission lines in the first overtone band of CO in the magnetic Herbig Ae star HD 101412. 
High-resolution CRIRES spectra reveal unusually sharp features which suggest the emission is formed in a thin disk centered 
at 1 AU with a width 0.32 AU or less. A wider disk will not fit the observations. Previous observations have reached similar 
conclusions, but the crispness of the new material brings the emitting region into sharp focus. 
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1. Introduction 

| The star HD 101412 (V1052 Cen) has been referred to in nu- 
merous papers as a Herbig Ae star. It has an infrared excess, 
and a disk (cf. Fedele, et al. 2008). 

HD 101412 is most unusual in having resolved, magnet- 
ically split spectral lines which reveal a surface field modu- 
, lus that varies between 2.5 to 3.5 kG (Hubrig, et al. 2010). 
Salyk, et al. (2011a) have surveyed molecular emission in a 
variety of young stellar objects. They found the emission to be 
much more subdued in Herbig Ae/Be stars than their cooler 
congeners, the T Tauri stars. This was true for HD 101412 as 
well, which was among the 25 Herbig Ae/Be stars they dis- 
cussed. One exception, however, was the molecule CO2, which 
had a very large flux in HD 101412; indeed, only one T Tauri 
star had a higher CO2 flux. 

In this letter, we report observations of lines in the first 
overtone (v=2 — > 0) band of carbon monoxide (CO) in HD 
101412. This feature has been the subject of previous studies 
(cf. Najita 1996, and references therein), leading to conclusions 
in general agreement with those presented here. 

Wheelwright, et al. (2010) studied this band in 7 upper 
main-sequence stars with masses from 6 to 43 M , signifi- 
cantly larger than that of HD 101412, ca. 2.5 M Q (Hubrig, et 
al. 2009). Individual lines were not resolved, but models were 
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* Based on observations obtained at the European Southern 
Observatory (ESO programme 087.C-0124(A)). 



constructed to fit the band head. The fits constrained the CO 
emission within limits of a few tenths of an AU (inner) to up to 
8 AU (outer). 

2. Observations and reductions 

We present observations made with the ESO CRIRES spec- 
trometer (Kaufl, et al. 2004) using a slit width of 0.2 arc sec- 
onds. The achieved resolution was determined from weak tel- 
luric lines, and is slightly above 90000. The signal to noise 
is difficult to obtain empirically because of strong and myriad 
weak telluric features. Local tests of the rms fluctuations in 1 
to 3 A regions ( 10 to 30 pixels) yield S/N between 100 and 
200, averaging 153. 

We discuss observations in two CRIRES regions, 2290- 
2301nm, and 2303-2313. 5 nm, made on 5 April 2011 at rota- 
tional phase 0.94 (Hubrig, et at. 201 1). The first of these regions 
contains the band head of 12 C 16 (henceforth, CO will mean 
this isotopologue). This head is located at 2292.8976 nm, the 
position of the rotational line R(51): (v" = 0, J" = 51 — » V — 
2, J' = 52). [Note: The double prime refers to the lower level.] 

Stellar (rest) wavelengths were determined using telluric 
and identified atomic stellar features, primarily C 1, Mg 1 and 
Si 1 lines. These lines were measured in other CRIRES regions, 
covering 1065-1105 and 1585-1550 nm. The transformation 
to rest stellar wavelengths was then applied to the CO regions. 

Telluric features were removed using hot standard stars 
with large rotational velocities, observed at comparable zenith 
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Fig. 1. The black spectrum shows CO emission lines 
R(20)(right) to R(27)(left) in HD 101412 after removal of tel- 
luric lines using the standard star, HR 4537, j Cen. (red). 



distances. For the CO features, this procedure introduced both 
noise, and artifacts; it was necessary to judge the position of 
the continuum in both the target and standard stars. The min- 
imum of the central M-feature is particularly sensitive to this 
subjective judgment (see FigJTJ. 
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Fig. 2. Observed CO R(20)-R(27) in black. The observations 
are those of Fig.[T](black). We subtracted unity, and scaled the 
remainder vertically to fit the observations. R-branch labels are 
written above the M-shaped profiles. Rest wavelengths of these 
lines correspond to the inner vertex of the M's. R(19) is at the 
extreme right, but only the violet edge of the observed profile 
is seen. Several sharp spikes may be attributed to imperfect re- 
moval of the telluric features that may be seen in the standard 
star spectrum of Fig. Q] Red shows calculated molecular fea- 
tures for a temperature of 2500K. 



3. The CO emission lines and band head 

3.1. Optically thick or thin? 

We assume an optically thin emitting gas as a reasonable ap- 
proximation. This is based on the sharpness of the profiles 
shown in Figs. \T\ and [2] Self-absorption would smooth the 
sharp edges of the M-shaped profiles, filling in the centers, 
and spilling intensity over at the sides (edges) of the profiles. 
Additionally, we obtain a good fit to the band head using this 
assumption (Fig. [3}. We can also make an estimate of the op- 
tical depth from the excess CO emission above the stellar con- 
tinuum, using stellar parameters from Cowley et al. (2010) and 
Hubrig et al. (2009). The method is similar to that discussed by 
Salyk et al. (201 lb), though we use a toroidal geometry for the 
emitting region rather than a slab. Plausible assumptions lead 
to an optical depth of 0.04 for R(20). Details may be found at 
http://www.astro.lsa.umich.edu/~cowley/oldindex.html 



3.2. An empirical fit 

All observations discussed here are in the R-branch of the 
first overtone band of the (X 2 £ + ) ground electronic state of 
CO. Identifications of molecular features, and temperature- 
dependent predictions are based on the data file coxx.asc 
(Kurucz 1993). Only 12 C ie O lines were considered. All oth- 
ers were assumed to be too weak to be relevant, or outside our 
limited observational regions. 

Emission spectra (Figs. [2] and [3) were calculated with the 
help of Kurucz's data. We assumed that emission flux profiles 
can be represented by gaussian functions with central inten- 
sity / oc gfA~ 3 exp —XulkT, where gf and^-,, are the oscillator 



strength and the excitation potential of the upper level of the 
CO transitions. We adopted T=2500K. 

In many cases (see Section [U, one can determine a rota- 
tional temperature from the relative intensities of the individual 
lines (cf. Brittain, Najita, and Carr 2009). However, the rela- 
tive intensities of R(20)-R(27) hardly change from T = 2000 
to 3000 K, the range we believe appropriate (see below). A 
first calculation, using a single Gaussian profile for each line 
placed the emissions in the middle of the observed M-shaped 
profiles. Carr (2005) discusses the canonical, disk-model M- 
shaped profile. We eventually adopted a set of 21 Gaussians 
for each line, with 0.0666 nm FWHM (0.04 nm 1/e-width). 
They were weighted appropriately to reproduce the observed 
profiles (red, Fig.0. The Gaussian at the center, therefore had 
the lowest weight. Other choices for the FWHM and number of 
Gaussians could give similar results. The fit is thus cosmetic, 
apart from the central wavelength matches, and the overall pro- 
file widths of the calculated features. These match very well the 
molecular wavelengths and the expected widths (see Section|4]i. 

The velocity half-widths of the Gaussians, for a mean 
wavelength of 2309 nm is 47.8 km s _1 . This value, uncertain 
by w 2%, is in agreement with the Keplerian velocity 47.1 
km s oused below in a disk-model fit to the profiles. 

The strongest lines in the region of Fig.|2]are R(20)-R(27). 
However, higher-order lines, with wavelengths increasing for 
R(n > 51) fill in some of the region between the stronger fea- 
tures. These R(n > 51) lines diminish in intensity for increasing 
values of n. 

Fig. [3] shows the region of the band head. Values of 
/" are written over the centers of the M-shaped profiles. 
Because of convergence to the head, vertical "sides" of the 
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Fig. 3. First overtone band head with partially cosmetic fit. The 
assumed splitting and relative intensities (within a given M- 
shaped line) of the multiple Gaussian profiles are the same as 
in Fig. [2] The overall distribution of intensities is somewhat 
better fit by assuming T = 2500 K rather than 2000 or 3000 K. 
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Fig. 4. Model fit (red) to the R(20) profile shown in Fig. |2] A 
thin, vertical line (green) marks the rest wavelength of R(20). 
The model is virtually a single ring, since r; n = 0.96 and r out 
1 .04 Al' 



profiles approach one another, coincide, and eventually cross. 
Additionally, close to the head, transitions in the R-branch lead- 
ing "away" (to the red) from the head (e.g. R(52)..R(56)) over- 
lap and are of comparable intensity to the lines with R-values 
below that of the head, R(51). 

The badly-blended head is reasonably well fit by the multi- 
Gaussian model. This indicates that the assumed profiles re- 
main a good approximation for rotational lines somewhat 
higher than R(39). However, the relative importance of these 
higher-series lines diminishes rapidly. The R(51) line is 3.3 
times weaker than R(31) at 2500 K. 

The overall shape of the head region is more sensitive to 
the rotational temperature than the region of Fig. [2] By-eye fits 
to the observed features were best for temperatures between 
2000 and 3000 K. The plot is for 2500 K, with an uncertainty 
of several hundred degrees. 

4. A model emission ring 

In this section we discuss a physical model to reproduce the M- 
shaped profiles. The code is far simpler than those used in some 
cited studies (e.g. Salyk 2011b). The simplification is possible 
because we assume (1) the emitting gas is optically thin, (2) 
in pure Keplerian rotation, and (3) unobscured by the star or a 
dust disk. 

The relevant code calculates the profile of a single emission 
line, and has not been adapted for multiple emissions, such as 
those shown in Figs.|2]and [3] 

We constructed a two-dimensional model of the emission 
region, consisting of a set of 81 concentric rings, each with 
288 emitting gas elements. They are all in Keplerian rotation 
around a 2.5 M star. The rings may be placed at arbitrary ra- 
dial distances from the star. We assume that each element emits 
a Gaussian profile with a FWHM of 0.0666 nm. 

A model that marginally fits the observations, has rings cen- 
tered at 1 AU, spaced by 0.004 AU. This gives a disk 0.32 AU 



wide. A wider disk, at 1 AU, gives a profile that is too broad 
and rounded at the peaks to fit the observations. Indeed a single 
ring would produce an entirely satisfactory fit. 

According to Fedele et al. (2008) the disk plane is 10° from 
the line of sight. Since cos(10°) = 0.985, we ignore the incli- 
nation, apart from the assumption that it must be enough for 
the emitting ring to be visible. A cosine factor may easily be 
incorporated, should future work reveal a larger inclination to 
the line of sight. Contributions from the 81 x 288 elements are 
Doppler shifted, added linearly, and scaled vertically to fit the 
observations. 

Fig. [4] shows a fit of the ring model to the rotational line 
R(20). The 81 rings are centered at 1 AU, and are separated by 
0.001 AU. 

Salyk et al. (2011b) stacked a number of their profiles in 
order to increase the signal to noise. We have attempted this in 
Fig. [5] The abscissa is rendered in velocity space. Plots have 
been scaled vertically for an optimum fit. The overall shape 
of the composite profile (black) is well fit by the model, apart 
from the low regions near the center of the 'M' (thick red). It 
is unclear how significant this is. We noted earlier the sensitiv- 
ity of this portion of the profile to subjective judgments of the 
continuum location. 

Additional emission is seen at the foot of the violet (or 
negative- V) side of the M-shaped profile. This is surely due 
to R-branch lines, R(7"), where J" > 5 1 . There is no provision 
in the model for this effect. 

Fig. [5] shows two additional curves which illustrate the ef- 
fect of an extended disk. Both curves assumed a disk with in- 
nermost rim at 0.4 AU, and outermost rim at 3.6 AU. The solid 
blue curve was made assuming each ring had the same total 
emissivity. The dashed green curve shows a profile where the 
emission is weighted by the negative 1.5 th power of the ra- 
dial distance. In this case, the innermost rings, with highest 
Keplerian velocities, have higher weights. This causes the ad- 
ditional emission at larger positive and negative velocities. 
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Fig. 5. Model fit to a composite of R(20) through R(27) of 
Fig. [2] The same parameters used for the fit of Fig. [4] (red) give 
excellent agreement for the composite. Two additional curves 
illustrate profiles that would result from a wider disk (see text). 



It is clear that such a disk cannot possibly account for 
sharply-edged rotational profiles. Quite elaborate radial inten- 
sity variations were considered by Acke, et al. (2005, cf. their 
Fig. 21), in connection with [O I] emission. Happily, such com- 
plexity is unnecessary in the present case. 



5. Discussion 

Van der Plas et al. (2008), and Fedele et al. (2008) have dis- 
cussed the protoplanetary disks around "three young interme- 
diate mass stars," including HD 101412, based on observations 
of forbidden oxygen, [O i] /16300. They discuss disk models, 
though not as localized as the one that fits our CO observations. 

Salyk, et al. (2011b) studied the fundamental (v=l -> 0) 
CO band in 31 T Tauri and Herbig Ae/Be stars, including a 
number with transitional disks. Rotational lines of the funda- 
mental CO band are more commonly seen in T Tauri stars 
than the first overtone, the subject of the present paper. Salyk 
et al. conclude their line fluxes are "consistent with emission 
from a single temperature ring width 0. 15 and 0.01 AU respec- 
tively...". This conclusion is very close to that of the present pa- 
per. However, the profiles illustrated by Salyk, et al. are much 
less sharply defined than those we observe in the first CO over- 
tone in HD 101412. The central portions of their profiles are 
mostly filled in, indicating emission from regions of a disk be- 
yond the innermost rim. Thus, they were also able to obtain fits 
(see their Table 4) with disk models extending over several and 
even tens of AU. The HD 101412 observations are incompati- 
ble with emission from such an extended disk. 

Perhaps the T Tauri star with CO (fundamental) profiles 
most closely resembling the first overtone lines of HD 101412 
is V836 Tauri (Najita, Crockett& Carr 2008). But note the slop- 
ing sides of their averaged profiles (their Figs. 4 and 5) in com- 
parison to the nearly vertical sides of an analogous profile in 
HD 101412 presented below (Section @). 



The existence of narrowly-confined emission regions is 
well established. The current observations may be the best- 
defined example found thus far. 
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